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Abstract. The matter distribution around the Perseus-Pisces 
(PP) supercluster is studied by comparing peculiar veloci- 
ties given by the Kinematics of the local universe (KLUN) 
galaxy sample to those predicted by Tolman-Bondi (TB) mod- 
els. To restrict the TB solutions we first solve the mass of 
the densest part of PP. This part is identified as a sphere at 
(l,b) = (140.2°, -22.0°), d « 50 h" 1 Mpc having a radius 
of 15/i _1 Mpc. This sphere surrounds the main part of the 
PP ridge and four most prominent clusters of the region. Us- 
ing virial-like mass estimators we calculate the cluster masses 
and obtain the upper and lower limits for the mass inside the 
15 ft,- 1 Mpc sphere: M PP =4 - 7/i _1 lO 15 M . This corre- 
sponds to a mass overdensity <5pp w 4, or ppp = 1-2 p cr . 
Mass to light ratios of the clusters are M/L = 200 - 600 
HMq/Lq, giving mass density ratio fin = 0.1 - 0.3, if the 
value of Mj L is assumed to be representative elsewhere in the 
universe. 

We estimate a radial density distribution around the PP core 
using two toy models and a smoothed density distribution ob- 
served for IRAS galaxies. The cosmological density parame- 
ters (fin, Oa) an d the PP mass are free parameters in the TB 
calculations. The KLUN velocities, obtained by Tully-Fisher 
relation and the normalized distance method, are adjusted by 
the Local Group (LG) infall velocity towards PP. Compari- 
son of the TB velocities to KLUN data points indicates that 
the infall velocity v m f < 100 km s . Allowing Mpp to vary 
within the limits given above we get constraints for the value 
of fio; fin = 0.2 - 0.4 are prefered to the more extreme values, 
fin = 0.1 or 1. A choice of either fiA = 1 — fio or do not 
cause any significant changes in the results. 

The validity of the TB model in complex environments is 
studied with an N-body simulation. There we see that the radi- 
ally averaged velocity fields around simulated clusters are com- 
patible with the corresponding TB velocities. This confirms the 
applicability of the TB model around large galaxy concentra- 
tions, providing that smoothed density and radially averaged 
velocity fields are used. 



Key words: Galaxies:clusters:general - Galaxiesxlusters: 
Pers-Pisc supercluster - Cosmology:theory - large-scale struc- 
ture of the Universe 



1. Introduction 

Tolman-Bondi (TB) model is the analytical solution of a spher- 
ically symmetric mass distribution in an expanding universe 
(Tolman 1934, Bondi 1947). It has been succesfully applied 
to the Virgo system by Teerikorpi et al. (1992) and, recently, 
by Ekholm et al. (1999) for describing the peculiar velocity 
field and estimating the total mass of the system, even though 
the galaxy distribution in Virgo deviates notably from spheri- 
cal. Motivation for the use of the TB model is that it is simple 
and analytical. It is a first step towards approximating an in- 



homogeneous universe. Bondi (1947) himself wrote that "the 



assumption of spherical symmetry supplies us with a model 
which lies between the completely homogeneous models of 
cosmology and the actual universe with its irregularities". 

TB results are valid with appropriate smoothing; the ob- 
served velocity field must be averaged over concentric shells 
around a galaxy concentration, when compared to the TB re- 
sults. Thus the TB model does not give the results in as much 
details as the various more commonly used gravitational insta- 
bility - reconstruction methods. On the other hand, it does not 
require any assumptions for the breakdown of the theory in the 
non-linear region. 

In order to test whether the TB model may be applied 
to an apparently still less favourable situation than Virgo, we 
study the Perseus-Pisces supercluster, where a TB-like devia- 
tion in the Hubble flow has already been observed (Theureau 
et al. 1998a). Our criterion for the success of the TB method is 



whether the mass inferred using TB is comparable to virial and 
luminosity x M/L estimate for the region, and whether the 
peculiar velocity field is adequately described by the model. 

The local d < 70h^ Mpcf] universe is dominated by two 
large mass concentrations, the Great Attractor and the Perseus- 
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Pisces (PP) complexes. According to Sigad et al. (1998) and 
Dekel et al. ( |1999| ) both regions have overdensities 5 = 1- 
1.5, where 5 = [p(x) — p]/p, p{x) is the local, and p the aver- 
age density. Both studies used a Gaussian smoothing of 12 h~ x 
Mpc and assumed f2o = 1. Thus the velocity field surround- 
ing PP should reveal a strong infall towards the supercluster's 
center of mass. For a galaxy at around half way from PP to us 
the infall velocity should be w 500 km s _1 . However, the in- 
fall signature is difficult to reveal because of large uncertainties 
in redshift independent distance measurements, e.g. those pro- 
vided by Tully-Fisher (TF) or Fundamental Plane (FP) meth- 
ods. Uncertainties are caused by statistical biases difficult to 
correct for in a heterogeneous environment. Additional inac- 
curacies are raised by the filamentary shape of the structure, 
consisting of an alignment of several subclusters, which com- 
plicates considerably any attempt to model the velocity field. 

1.1. Previous studies of Perseus-Pisces 

In the past ten years PP has been investigated by numerous 
authors using various distance criteria and statistical methods. 
The region is now well covered by great amount of data. From 
a large reshift sample of 3311 galaxies in the region bounded 
by 22 h < R.A. < 4'\ +20° < Dec. < +50°, Wegner et al. 



(1993) described the Perseus-Pisces complex as follows: The 



supercluster consists of a continuous arrangement of high den- 
sity clusters and groups, among which the Pisces, A262, and 
Perseus (A426) clusters are the most prominent features. Its 
main ridge extends at least 50 hr 1 Mpc from Pegasus eastward 
to the Perseus cluster, with a width in the plane of the sky of 
about 5-10 hr 1 Mpc, and a typical depth in redshift of 250- 
500 km s _1 . All the objects in the ridge are lying roughly at 
the same distance from us, at ~ 5000 km s _1 . 



First kinematical studies (Willick 1990, 1991, Courteau et 



al. 1993) based on an r-band TF sample of 355 galaxies, sug- 



gested that the whole structure takes part on a large scale flow 



defined by the Great Attractor (see e.g. Lynden-Bell et al. 1988, 
Dekel 1994), with a mean velocity of ~ 350 km s _1 towards 



based on Monte-Carlo simulations, in order to control biases 
and gain a better sampling of the sky. Their studies revealed a 
clear infall motion towards PP, while the GA was found less 
dominant than originally believed. 

The differences in the results stated above may be caused 
by uncertainties in the reconstruction method, in the correc- 
tions for the selection effects, and/or in the distance indicator 
calibration. For example, an error in the TF zero point would 
cause a systematic error that rese mbles a Hubble flow in the 
peculiar velocity field. Dekel et al. ( 1999 ) also note that the re- 
sults in the PP region are more sensitive to noise than the other 
regions. 

1.2. What is new in our approach to PP? 

In a previous study using the KLUN galaxies, a full infall pat- 
tern, both in the front and the backside of the PP supercluster 
was shown, for the first time. A new Malmquist bias correc- 
tion method for the TF distance data was used (Theureau et 



al. 1998a). In the present paper we investigate the dynamics of 



galaxies in the PP surroundings in more detail. We study the 
spatial distribution and peculiar velocities of the galaxies using 
the Malmquist corrected TF distances. The observed velocity 
field is averaged over spherical shells around the PP center and 
compared to the results given by TB calculations. 

The outline of the paper is as follows: In Sect. [| we 
briefly present the data and summarize the normalized distance 
method (NDM) that has been used for deriving unbiased aver- 
age distances and peculiar velocities. In Sect. || we calculate the 
virial masses of the most prominent clusters in the PP region 
and estimate the mass of the whole central region of PP from 
these masses. Section [| includes a description of the Tolman- 
Bondi model and its application to PP. The model is a lso tested 
on an N-body simulation of galaxy clusters (Sect. 4.2). In Sect. 

we compare the TB model to the infall pattern observed with 
the KLUN data. Finally, the main results are summarized in 
Sect. H 



the Local Group (LG). Han & Mould ( |1992[ ) confirmed this 2 . Towards unbiased peculiar velocities 



result independently using a whole sky TF sample of 21 clus- 
ters. Hudson et al. ( 1997 ) combined their inverse FP relation 
measurements of seven clusters in the PP region with data of 
nine other clusters from the literature and concluded that the 
whole region has a mean CMB frame bulk motion of 420 ± 280 
km s _1 towards galactic coordinates (l,b) = (262°,— 25°). 
POTENT reconstruction method of the densit y and velocity 
field for Mark III galaxies (Dekel et al. 1999 ) gave different 
results. There it was found that the center of PP is roughly at 
rest in the LG frame and the regions between the LG and PP 
are infalling towards PP. 

Most studies of the velocity field of galaxies suggest that 
the kinematics of the local universe can be described by a 
model in which galaxies are infalling to two mass concen- 
trations , one in PP and on e in H ydra-Centaurus (e.g. Han & 

1999t) . However, Freudling et al. 



Mould 1992, Dekel et al. 



The following notation is appointed for different distance and 
velocity variables: 

- d for distances from the LG 

- r for distances from the PP center 

- capital R for fixed distances: i?pp(= 50/i _1 Mpc) is the 
distance between the LG and the PP center, R s (= 15/i _1 
Mpc) is the radius of the main region of the PP concentra- 
tion 

- small v for peculiar velocities, v r is the peculiar velocity 
with respect to the PP center, Vd is the peculiar velocity 
with respect to the LG 

- capital V for observed (Hubble flow + peculiar) velocities, 
corrected for the infall towards Virgo as described in Sect. 
2T 



(|1995j) and da Costa et al. ( |1996| ) used a method of analysis, Exceptions are 
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- quantities with subindex 25: D25 is the standard notation 
for the isophotal B = 25 mag diameter of a galaxy and 
i?25 is the corresponding axis ratio 

- V m is the maximum of the rotational velocity of a galaxy 



Direct (B-magnitude) TF relation 



For the Hubble constant we 



km s 



Mpc~ 



or indicate it 



assume Hq = 60 
by the term h = 



i?o/(100 km s 1 Mpc 1 ). The value is needed in calculating 
the radial part of the peculiar velocities, 

Vd,obs = V - H d, (1) 

where the subscript "obs" emphasizes that only this component 
of the peculiar velocity can be directly observed. 

The main difficulties in velocity field studies are caused by 
the large uncertainties in distance determination and by the re- 
sulting biases. Attempts to overcome these problems require 
a few things. Firstly, one must know thoroughly the sample, 
how it was collected, and to which extent it can be assumed 
complete. Secondly, systematic effects in the distance criterion 
have to be modeled, so that one can derive for a given appli- 
cation an unbiased set of distances. In particular, biases will be 
different if one estimates peculiar velocities at a given derived 
distance, at a given true distance, or at a given observed radial 
velocity. Thirdly, the sample has to be deep enough and should 
cover a sufficiently large solid angle in the sky for the average 
relations, needed for the bias correction, to be free from local 
perturbations in the density or the velocity field. This last re- 
quirement is similar as needed for a satisfactory determination 
of the Hubble constant. 

2.1. The data 

We use the KLUN (Kinematics of the local universe) sam- 
ple consisting of 6600 spiral galaxies having measurements of 
isophotal diameter D25, HI line width, radial velocity, and B- 
magnitude. This sample has been fully des cribed in previous 
papers of this series, e.g. by Theureau et al. (1997), so we only 
summarize here the main characteristics. 

KLUN is selected according to apparent diame ter, and is 
complete down to D25 = L'6 (Paturel et al. 1994 ), covering 
the type range Sa-Sdm (T=l-8). The data are extracted from 
LEDA (Lyon-Meudon Extragalactic Database, Paturel et al. 



1997) and complemented with optical redshifts measured with 
ESO and OHP telescopes, and HI spectra obtained with Nangay 
and Parkes radio telescopes (Bottinelli et al. 1992, 1993, di 



Nella et al. |1996| , Theureau et al. |1998b[ ). According to LEDA 
precepts, all astrophysical parameters have been homogenized 
and reduced to a standard and common system (Paturel et al. 



1997). Isophotal D25 diameters and apparent B -magnitudes are 
corrected for galactic extinction according to Fouque & Paturel 
(1985), and for inclination, i.e. the opacity effect, as explained 



in Bottinelli et al. (1995). HI line widths, reduced to the stan- 



dard levels of 20% and 50%, are corrected for internal velocity 



dispersion according to Tully & Fouque (1985). Heliocentric 
radial velocities are change d to th e rest frame of the LG cen- 
troid following Yahil et al. ( 1977 ) and corrected for a Virgo- 
centric component, assuming an infall of the LG towards the 



Unbiased plateau 













Fig. 1. The logiJo vs. d n diagram of KLUN galaxies. The 
line, indicating the average (log_ffo)(^n), is (approximately) 
horizontal in the unbiased region at small d n . Malmquist bias 
makes the curve to rise at larger d n . 



Virgo cluster vq = 200 km s (Theureau et al. 1997, Jerjen 
& Tammann 1993) and a radial velocity of Virgo Vb,vir = 980 



bles 1976) 



km s 1 (Mould et al. 1980), within a linear infall model (Pee- 



In the TF relation studies we need to make some restrictions 
to the data. Galaxies close to the Galactic plane (|&| < 15°) are 
excluded due to the large uncertainties in the Galactic extinc- 
tion correction. Face-on galaxies, log R25 < 0.07, are excluded 
to prevent large errors on calculation of the rotational velocity 
parameter p = log V m , where V m is the maximal rotational ve- 
locity of a galaxy, measured from the HI line width. 

Absolute TF distances are obtained using the calibration 



given by Theureau (1998). TF scatter was there reduced by 
30% by taking into account the dependence of the TF zero 
point on mean surface brightnesses of the galaxies. The NDM 
first order Malmquist corrections are calculated on the basis of 



the whole sky sample, according to Theureau et al. (1998a). 
The sample with thus corrected TF distances, restricted to the 
strictly diameter complete part, contains finally 2800 spiral 
galaxies. 

2.2. The normalized distance method (NDM) 

The concept of normalized distances was introduced in a theo- 
retical discussion by Teerikorpi (1984). The method was first 
applied by Bottinelli et al. ( 1986| ) and improved in a series 
of papers, most recently by Theureau et al. (1997). A theo- 



retical description in terms of probability densities was given 



by Theureau et al. (1998a), where the authors showed that the 
method can be extended to make Malmquist bias corrections to 
galaxy distances even beyond the so-called unbiased plateau. 
The unbiased plateau has been profitably used for determina- 
tion of the unbiased value of Hq and has the advantage that one 
does not need model-dependent bias corrections. 

In order to understand how the NDM works, let us define 
classes of galaxies as follows: objects belong to a given class if 
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they have identical values of log rotational velocity p, Galactic 
extinction correction a g , inclination R25, morphological type 
T and mean surface brightness S. The main idea of the NDM 
is that for a given magnitude or diameter limited sample, and 
for a given TF scatter the behaviour of the Malmquist bias^] as 
a function of true distance is similar for all galaxies belonging 
to the same class. The class is assumed to have a Gaussian dis- 
tribution of absolute magnitudes. Then, for each class, there is 
a unique curve of evolution of the bias, which is conveniently 
expressed as the curve of observed average value of the Hub- 
ble parameter (log Ho) = (log V c — logd^p) vs. distance (d), 
where c£tf is the distance given by TF relation, 



log^TF = otfP + ^tf(S) - log-D 2 5 - log 



108' 



otf is the TF slope, &tf(£) the TF zero point, dependent on 
the surface brightness £ of the galaxy, and is a constant for 
converting the units (linear TF-diameter in kpc and D25 in O'.l 
to c?tf in Mpc). 

Essentially, the normalized distance technique means su- 
perimposing all the bias curves (log Ho) (d) of the different 
classes of objects by using an appropriate scaling on the dis- 
tance: 



dn = d- f(p, a g ,R 25 ,T,Y,) 



(3) 



As a result, the amplitude of the bias can be expressed merely 
as a function of the normalized distance (d n ), the apparent limit 
of the sample (mu m or Z?25,iim), and the TF scatter (cttf)- Fig- 
ure IT] shows the (log Ho) vs. d n diagram of KLUN galaxies in 
the case of the B magnitude TF relation. The unbiased part of 
the sample is identified as a plateau of (log Hq) = constan- 
tat, at low values of d n . Note, however, that for calculation of 
d n , one needs the true distance d, usually derived as the so- 
called kinematical distance, e.g. in the case of the Hubble law 
d = dkin = V/H Q . 

The average bias correction is valid only if the following 
assumptions are fulfilled: 

- HI : the sample is complete and strictly diameter (magni- 
tude) limited 

- H2 : TF residuals follow a Gaussian distribution, and the 
TF scatter does not depend on the value of p 

- H3 : the selection in apparent diameter (magnitude) is in- 
dependent of the selections in redshift and p 

- H4 : the radial velocity field integrated over all the direc- 
tions vanishes, i.e. except the Hubble flow, there is no co- 
herent flow at the scale of the whole sample 

Assumption HI is obtained by applying a strict cut-off 
at the completeness limit of the sample. H2 is in agreement 
with what is observed in the case of the B-band TF relation 
(Theureau et al. 1997 1. H3 might be difficult to prove, if obser- 



vational errors on apparent magnitudes or diameters are corre- 
lated with the errors on p : both are linked to the inclination of 



In fact the Malmquist bias of the 2nd kind, as termed by Teerikorpi 



the galaxy through the internal extinction correction and the de- 
projection of the observed HI line width. However, this effect 
is insignificant in the case of the diameter TF relation. 

As for assumption H4, a list of recent measurements of 



large-scale bulk motions is given by Strauss (2000). The list 



is somewhat divergent; roughly half of the surveys find no bulk 
flow at scales V > 6000 km s -1 , the rest claim bulk flows of 
300-600 km s _1 . We will return to the question of large-scale 
flow within the KLUN sample in the future, for now we assume 
that assumption H4 is true for our data. 

Finally, the corrected distance is given by 



(2) log d c = log d TF 



2 &tf exp 



(- 



(qj-log_D 2 5,lim)' 



1 - erf 



/ M — log £>25,lim A 

\ \/2cttf / 



, (4) 



(Theureau et al. 1998a 1, where the (log) normalized distance is 
marked with u>, 

ui = log d n = \ogd + dTF (2.7 - p) + 0.094a g 
-Clogi? 25 -&T F (£). 

Here d — V/Hq and parameter C defines the correction for in- 
clination (Bottinelli et. al, 1995). For a single galaxy class and 



usual kinematical distance, an analogous formula was derived 



by Teerikorpi (1975). It has also appeared in the "triple-value 



correction method" of Sandage (1994 ) 



We emphasize that this correction is statistical in nature. 
When applied to an individual galaxy, it gives an accurate av- 
erage correction only if the velocity field model (providing c?kin 
and d) is correct. When the velocity field is distorted, say, in the 
direction of a supercluster, it still gives a first order correction, 
but a com plete correction would require an iterative procedure 
(see Sect. 53). Even the first order correction essentially im- 
proves thederived velocity field, as was shown by Theureau 
et al. (1998a I for the PP supercluster, but it is notable that the 
peculiar velocities will be somewhat suppressed if only a first 
order correction is done. 

3. The virial mass of Perseus-Pisces 

Figure || shows a map of the PP region. Circles stand for galax- 
ies from LEDA, having radial velocities between 3000 and 
7000 km s _1 . This limitation allows us to look at the main 
part of the PP concentration, lying atVw 5000 km s _1 , prac- 
tically free from any contamination from fore/background ob- 
jects. LEDA is regularly updated for all published data and thus 
contains all the recent PP-surveys. 

Symbols other than circles are KLUN galaxies. +:s are ob- 
jects within the velocity interval stated above, *:s have V > 
7000 km s" 1 and x:s V < 3000 km s _1 . Only the KLUN 
galaxies satisfying the completeness limit, D25 > 1 '6, and 



(1997) 



other criteria demanded by our TF-method (see Sect. 2.2), are 
plotted. In total, Fig. || contains 2016 LEDA galaxies and 266 
objects from KLUN. 

In Fig. ||the Hubble diagram of the KLUN PP region galax- 
ies is shown. The line corresponds to the Hubble law with 
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Table 1. Galactic coordinates, numbers of galaxies, mean ve- 
locities, virial radii in degrees and in h^ 1 Mpc, upper and lower 
estimates of masses in /i _1 10 14 M Q , and mass to light ratios in 
HMq/Lq, of the four largest PP clusters in radial velocity in- 
terval 3000-7000 km s" 1 . Clusters are 1) Perseus, 2) A262, 3) 
0122+3305, and 4) Pisces. 
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Fig. 3. The Hubble diagram of KLUN galaxies in the PP region. 
d c is the NDM corrected TF distance. 



Hq = 60 km s _1 Mpc -1 . The infall towards PP is seen as 
the surplus of galaxies above the Hubble line at small distances 
and below the line at large distances. 

We seek a rough estimate for the amount of matter in the 
PP region. First we determine the masses of the densest clus- 
ters in the region. We define the densest consentrations as the 
ones having the projected galaxy overdensity > 10. Here 
<52d(x) = (p2b{x) — P2d)/P2D, the two dimensional sky den- 
sities p2u(x) are calculated in 1 degree 2 regions, and p2D is 
the mean value of the whole region shown in Fig. ^. The lim- 
iting factor 10 has no special physical meaning. It was chosen 
to give a reasonable number of dense galaxy concentrations. 
There are four clusters satisfying this criterium. The charac- 
teristics of these clusters are listed in Table [j]. Masses are esti- 
mated using the virial theorem and the alternative methods pre- 



sented by Heisler et al. ( 1985 ), see Table || below. The effects 
of possible outsiders and interlopers are tested by using various 
radii and redshift intervals in defining cluster members. Table 
U shows some of the results. For the mass evaluations in Ta- 
ble we adopt ^-interval of 3000-7000 km s - 1 , and use virial 
radii of 2° for Perseus and A262, and 1.5° for 0122+3305 and 
Pisces, These radii are drawn in Fig. ||. 

For the M/L ratio derivations in Tables [l] & [| we calcu- 
late the total luminosities of the clusters from LEDA B-band 
brightnesses and estimate selection effects using the Schechter 
function 



cj>(M) OC I00-MM>-M)( a + 1) exp [_ 10 0.4(M*-M)] 



(5) 



expressed here in terms of absolute magnitude Ai. The char- 
acteristic magnitude, M*, and slope, a, of the function are 
found by fitting Eq. |^ to the observed magnitude distribu- 
tion. Since the Perseus cluster shows clearly different fea- 
tures from the other three clusters, the fitting for it is done 
separately. Obtained parameters are a — ~l-46lg'2g anc l 
M* - 5 log h = -19.6^3 for Perseus, and a = -1.51^;^ 
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Table 2. Virial, projected, median and average masses by meth- 
ods of Heisler et al. (1985), and the corresponding M/L ratios 
of the four PP clusters. Different radii 9 V are used for each. The 
figures indicate that virial and median mass estimators seem to 
be less affected by outsiders (dM / dd v small) than projected 
and average masses. Projected mass is also very sensitive on 
the coordinates chosen for the cluster center. This table is cal- 
culated for galaxies having 3000 km s^ 1 < V < 7000 km s -1 . 
Changing the V interval affected all the mass estimators sim- 
ilarly. The limits for V and the cluster virial radii are chosen 
partly by eye, simply by using limitations defining a compact 
galaxy concentration, and partly by requiring certain level of 
stability for the obtained values of mass. The values adopted 
in Table [j] are the lowest and the highest values from the rows 
V = 2.0° for Perseus and A262, and V = 1.5° for 0122+3305 
and Pisces. Units: [M] = fo- 1 10 14 M©, [M/L] = hM /L Q 
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Fig. 4. Figure a. shows the observed absolute magnitude his- 
tograms for the Perseus galaxies (gray bins) and the three 
other clusters combined (hatched bins). Two Schechter func- 
tion curves are drawn on both of these histograms. Lower lines 
are fitted to the histograms at M. — 5 log h < —18, upper lines 
are calculated by assuming completeness of 80% for Perseus 
and 85% for the other clusters, following the latitude depen- 



dent completeness estimate by Marinoni et al. (1999). Figure 
b. shows the corresponding total luminosities in 0.5 magnitude 
bins. The lines are the values predicted by the upper Schechter 
curves of Fig. a. By estimating the ratio of the gray/hatched ar- 
eas to total areas under the Schechter curves in Fig. b. we can 
calculate the amount of luminosity of non-detected galaxies. 
This is 42+2% for Perseus and 40+J% of the total luminosity 
for the others. 



-20.41^ for the others (Fig. 



The 



and M * - 5 log ft, 
errors are bootstrap 1 a confidence limits. 

Marinoni et al. ( 1999| ) observed for galaxies in large groups 
a = -1.28 ± 0.18 and M* - 5 log h = -20.38 ± 0.31. The 
difference between our value for Perseus M. * and the one given 
by Marinoni et al. is probably due to environmental differences 
between Perseus and other clusters and/or uncertainties in the 
brightness measurements caused by the Galactic absorption. 
Comparison of the total observed luminosities of the clusters 
to the values obtained by integrating Eq. [| indicates that the 
contribution of non-detected galaxies is V n d = 42 + J% for 
Perseus and V n d — 40+3% f° r tne om ers, with 1 o bootstrap 
errors. The limits for the M/L ratios are then calculated with 
the upper and lower mass values from Table [j] and by estimat- 



ing the true total luminosities to be 1/(1 — V a d)> i- e - 1.64-1.80 
(Perseus) and 1.60-1.95 (others) x the observed values. Re- 
sulting M/L ratios are all in the range of 210-550 JiMq/Lq, 
in good agreement with other recent estimations of M/L ratios 



of rich clusters (e.g. Bahcall et al. 1995 Carlberg et al. 1997). 

Adding up the masses of the four clusters in Table |lj we get 
~ 2-3 /i~ 1 10 15 Mq. The center of mass of PP, using these clus- 
ters, is at V c = 4965 km s _1 , (fc,M = (140.2°, -22.0°). 
Distances of the clusters from this point are 11.5, 3.9, 9.9 and 
12.5 h~ x Mpc. A sphere at (lc,bc), with a radius of 15/i _1 
Mpc surrounds all these main concentrations. In addition to the 
four clusters, this sphere contains 513 other galaxies. If we as- 
sume that the M/L ratio does not grow in scal es larger tha n 
the rich clusters, as suggested by Bahcall et al. ( 1995 , 2000), 
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we can then estimate Mpp w 4-7 /i _1 1O 15 M0, where Mpp is 
the mass inside the 15 hr 1 Mpc sphere. The selection effects 
are assumed to be the same for the galaxies in the four clusters 
and for the other galaxies in the 15/i _1 Mpc sphere. 

However, there is some evidence that M/L increases on 
scales larger than clusters. For example Small et al. (1998) 
found the mass to light ratio to be 564 IiMq/Lq for the 
Corona Borealis supercluster and 726 h M Q /L Q for a back- 
ground supercluster. Small et al. used the virial and projected 
mass estimators for the whole bound supercluster. The Tolman- 
Bondi calculations (next Sect.) show that the PP is also bound, 
so we may apply the virial estimators on the PP also. We get 
Mpp = 2.6, 4.2, 2.8, and 3.3 r^O^Mg for the 836 galax- 
ies in the central 15/i _1 Mpc sphere in PP, using the virial, 
projected, medium, and average mass estimators, respectively. 
The corresponding sum of the galaxy luminosities gives 4.4 
Ji~ 2 10 12 Lq. Taking from above the more conservative limits 
for the number of non-detected galaxies, this number is then 
increased by a factor 1.60-1.95, giving the total luminosity 
L P p = 7.0-8.5 h~ 2 lQ 12 L Q . Finally, the mass to light ratio, 
using these values, is 330-600 HMq/Lq. However, since this 
mass for the whole PP core is about equal to the sum of the 
four clusters derived above, we doubt the accuracy of the virial 
and other estimators in this case, and adopt the values 4-7 
/^lO^AfgforMpp. 

These upper and lower limits for Afpp give a mean density 
Ppp » 2^1 /i 2 10 -26 kg m~ 3 « 1-2 p cr , where p CT is the criti- 
cal density of the universe. In terms of mass overdensity this is 
<5pp ps 4, assuming Slo ~ 0.2-0.4. The lower value of f2n here 
corresponds naturally to the lower limit of ppp, and high fig to 
the higher ppp. The result is in agreement with the value 6 « 1 
(O = 1) by Sigad et al. (|1998|) and Dekel et al. fll999[), 



4. Tolman-Bondi method 

4.1. Theoretical background 

Tolman-Bondi (TB) equations determine the kinematics of a 
spherically symmetric mass concentration in an expanding uni- 
verse (see e.g. Olson & Silk 1979, and Ekholm 1996 for main 
features of TB, and Teerikorpi et al. 1992 Ekholm & Teeriko- 
rpi 1994, and Ekholm et al. 1999, 200C for some recent appli- 



cations). In practice we study the behaviour of a spherical layer 
of matter at distance r from the center of symmetry. Each of 
these shells have two dynamical constants. In Newtonian terms 
we can think of these constants as the mass inside the radius of 
the shell, 



M(r) 



4np(r')r' 2 dr' 



o 



and the energy per mass of the shell, 

, , 1 fdrV GM(r) Ac 2 r 2 
E{r) ~ ' 



2 V dt 



(i 



(6) 



(7) 



In view of the recent interest in Friedmann models with non- 
zero cosmological constant A, we have added its contribution. 



The energy term determines whether the mass shells will even- 
tually escape (E(r) > 0) or collapse to the center {E(r) < 0). 
The usual assumption that the mass shells do not overtake each 
other, i.e. there are no shell crossings, is followed here. 

We attempt to solve the present day velocities of the mass 
shells, v r (r), for the spherically symmetric mass distribution, 
with the initial condition that at some point in the past these 
mass shells have all been very close together. That is, we as- 
sume the big bang to have happened at t — 0. This moment is 
acquired from the underlying Friedmann model, by calculating 
the age of the Universe To = To(Ho, S7o, fi/v) for the set of 
cosmological parameters Ho, fto = po/ Per, and £Ia = Pa/ Per, 
where p CI = p is the mass density, and pA = 4^q- Then 
we find the value of the energy term E(r) for each present day 
radius ro, by demanding r(To, E) = ro and r(0, E) = 0, and 
integrating backwards in time. The integration is performed nu- 
merically for both r and v r with 



r = j2^m + ^ + E 



(8) 



(9) 



ensuring the conservation of energy, and 

GM(r) Ac 2 r 

This numerical approach is used because for A / model 
a simp le parameterized solution given in e.g. Olson & Silk 
(1979) does not apply. For A = 0, we checked the numerical 
calculations by the parameterized formulae. 

For t — ► 0, where the velocities grow considerably, the nu- 
merical methods confront some difficulties. Even with adapted 
stepsize algorithms we can not avoid all the problems. Using 
a fourth order Runge-Kutta method, we obtain inconsistent re- 
sults, r 7^ v r , when t < ^j. Fortunately our main concern is 
to calculate the dynamics at larger t, where the method proved 
to be robust. A comparison to the parameterized TB solution 
showed that the problems at small t did not have any effect on 
our resulting values of E(r). With the E(r) thus obtained, the 
present day velocities v r (ro)\E=E{r ) can be calculated. 

4.2. Testing the method on an N-body simulation 

The TB model is a rough approximation of the dynamics of 
galaxy clustering. The reader may suspect that the assumption 
of spherical matter distribution around a galaxy cluster is per- 
ilous. In reality galaxies seem to accumulate into filamentary 
and wall-like structures. However, an averaged galaxy density 
field may be well described with a spherical model; the IRAS 
map, using a 12 hr l Mpc smoothing, shows a roughly spheri- 



cal density field around PP (Sect. \3). And the distribution of 



dark matter may be even more spherical. 

How about the peculiar velocities? The many-body gravita- 
tional interplay determining the dynamics of galaxy clustering 
is a complicated issue, can it really be described by smooth ra- 
dial TB velocities? That may seem doubtful, but, again, we rely 
on the power of averaging. The TB results are compared only 
to the average radial velocities around the center of mass. 



8 



Hanski, M.O. et al.: Perseus-Pisces and Tolman-Bondi 













100 - 














7 






















50- 


















Y 






X ^' b0 

v^100 






: nop 





Fig. 5. The (141 h^ 1 Mpc) 3 cube of GIF N-body simulation 
galaxies. The simulation started with 256 3 dark halos, produc- 
ing the 15000 galaxies seen in this figure. The four circles in- 
dicate the clusters for which we studied the applicability of the 
TB calculations. 



The correspondence between the TB solutions and the av- 
eraged observational quantities is here studied with a cosmo- 
logical N-body simulation. In a simulated galaxy cluster we 
have the advantage of knowing the true distribution of mat- 
ter and the true velocities, not distorted by selection effects or 
observational uncertainties. And even better, we can study the 
evolution of the cluster and compare the matter and velocity 
distribution of the galaxies to those predicted by the TB model 
at different epochs. 

In what follows we use the N-body simulation 
data provided by the GIF group, a German-Israeli col- 
laboration working with the Virgo consortium project 



(http : / / star- 



.dur.ac. uk/ ^frazerp/virgo/) 



The GIF team investigates the clustering of gal a xies in 
a hierarchical u nivers e (Kauffmann et al. 1999a 1999b , 
Diaferio et al. 1999). They simulate structure formation 
in a periodical (141 hr 1 Mpc) 3 cube using a cold dark 
matter, flat curvature, non-zero cosmological constant 
model with il = 0.3, Q A = 0.7, and h = 0.7. The 
simulation starts with 256 3 dark matter halos with masses 
1.4/i _1 10 10 Mq. The galaxies are formed at concentrations 
of halos, and the effects of gas cooling, star formation, and 
supernova feedback are taken into account. The data files at 







Fig. 6. The four GIF galaxy clusters extracted for further study. 
Figures show XY-planes from the previous figure, centered to 
the cluster centers. The width in Z direction of each plot is 30 
Mpc. Circles with radii R s are drawn. It is notable that these 
complexes are not spherically symmetric, resembling the situ- 
ation with PP. 

Table 3. Radii R s , and masses and average overdensities inside 
i? s , of the four simulated clusters in Fig. |6| The units are [Rs] = 
h' 1 Mpc, [M] = h- 1 1O 15 M . 





R s 


M 


{S)n, 


1 


10 


4.5 


8.0 


2 


14 


7.2 


4.3 


3 


8 


3.0 


10.3 


4 


10 


4.4 


7.9 



Figure |] shows the galaxies in the whole (141 h^ 1 Mpc) 3 
cube at z = 0. We study the evolution of four large clusters in 
the simulation. These clusters are circled in Fig. || and shown 
individually in Fig. |[ For each of these clusters we approxi- 
mate the central radius R s by eye, and calculate the masses and 
the radially averaged S(r), see Table |[ These values are similar 
to what we observe in the PP supercluster. 

The overdensities are calculated from the halo number dis- 
tribution, with a 3/i _1 Mpc smoothing. At this smoothing 
length the matter distribution in the simulation is unbiased ac- 
cording to Kauffmann et al. (1999a). Then we determine the 



center of mass and calculate the radially averaged S(r) around 
http : / /www . mpa-garching . mpg . de /Virgo/ data_do6^ri(rf(lh£hffaimil{ ;lusters. For the TB calculations we use either 



have the positions, velocities, and luminosities of galaxies and 
distribution and masses of halos at six redshifts, z =0.00, 0.20, 
0.42, 0.62, 0.82, and 1.05, corresponding to epochs t/To = 
1.00, 0.97, 0.79, 0.67, 0.57, and 0.49. There are two sets of 
files, one for the galaxies (luminous matter) and one for the 
halos (luminous + dark matter). 



this observed S(r) (at z — 0) or a "step function" toy model 
where 5(r < R a ) = (S) Rs , and S(r > R s ) = 0. Figure ^ 
shows the evolution of S(r) in the simulation and by means 
of the TB model. The toy model is very instructive in showing 
how the density enhancements grow in amplitude and converge 
in comoving radius with time. 
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Fig. 7. The density fields S(r) around the four simulated clusters at the six epochs between the present and z = 1.05. Points show 
the averaged 5(r) of the simulation, calculated using a 3/i _1 Mpc smoothing. Error bars are the 1 a standard deviations. Gray 
line is the observed density field at z — evolved backwards in time via the TB model. Thick black line is the TB evolved "step 
function" toy model. R s are the adopted cluster radii and S s are the average overdensities within these radii at z = 0. 



As expected, in the simulation the structures form at z > 1 
and there is not much evolution in 6 (r) from z = 1 .05 to z = 0. 
Contrarily, in the TB evolved clusters the overdensities are low 
at z = 1.05 but grow quite steadily to the present day values. At 
larger r the S(r) given by TB model follow the corresponding 
values of the simulation quite well at each epoch. 

Figure || shows the comparison of the observed averaged 
peculiar velocities in the simulated clusters to the TB results. A 
similar time dependent factor as in the evolution of the density 
field is seen here. At earlier epochs the peculiar velocities given 
by the TB model are somewhat smaller than the v r (r) in the 
simulation. At z = the TB values are higher, although the 
differences are not at a significant level. 

Our conclusion is that concerning the evolution of density 
and velocity fields around clusters there is a marginal difference 
between the simulations and the TB calculations. The overden- 
sities 5(r) are quite well approximated at large r at each epoch, 
and the present day peculiar velocities v r (r) are sufficiently 
well derived, that is, with < 1 a deviation, using the observed 
present day S(r). 



4.3. TB method applied to Perseus-Pisces 

Now, how to use TB method for Perseus-Pisces where the 
distribution of luminuous matter apparently is not spherically 
symmetric? Certain restrictions and assumptions must be ap- 
plied. Firstly, we isolate a sphere containing the main part of 
the mass, as was done in Sect. |[ Then we assume that any sub- 
structures or dynamical features inside the sphere do not affect 
a galaxy outside it; the sphere can effectively be replaced by a 
point mass located at the center. Outside the sphere, we assume 
the density to decrease to the average density of the universe 
following spherical symmetry: 

8{r,6,(f>) = 5{r) if r > R s 
S(r) -> for r > R s 

(R s = radius of the sphere). After these assumptions we can 
apply TB equations for galaxies close to, but not inside the 
sphere. 

The assumptions may seem doubtful. Is it not oversimpli- 
fied to approximate the filamentary structure seen in the PP 
region with a spherical model? The filaments must produce 
observable non-radial motions. The main justification for our 
method lies in the power of averaging. Small scale non-radial 
peculiar velocities are smoothed away in a similar manner as 
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Fig. 8. The evolution of the velocities v r (r) around the four clusters. Crosses and error bars are the averaged v r (r) and their 
standard deviations observed in the simulation. Gray lines correspond to the observed density distribution and black lines to the 
toy model density distributions, evolved within the TB models. 



a smoothed density map blurs away small scale substructure. 



Compare the smoothed IRAS map by Sigad et al. (1998) to the 
galaxy map in Fig. ^| The IRAS map shows a spherical density 
distribution centered in the PP, quite unlike the first glance on 
the individual gala xy distribution suggests. N-body simulations 
indicate (Sect. \2) that the same thing happens for the velocity 
fields, justifying the usage of the TB models. In any case, over 
90 percent of the matter is dark and may be more symmetrically 
distributed. The final criterion is how well the model describes 
the infall pattern and predicts the total mass. 

For the density distribution outside the i? s -sphere we con- 
sider three models: two toy models, a step function 



Si(r) 



Spp if r < R s , 
if r > IL, 



a declining profile 

5 2 (r)=Ar-\ 



(10) 



(11) 



where A is chosen so that M(R S ) = Mpp, and a rough approx- 
imation of the smoothed density field of IRAS galaxies around 
PP, taken from Sigad et al. ( 1998[ ), where the middle panel of 
their Fig. 6 shows the supergalactic plane with LG in the cen- 
ter and PP at (X,Y) w (50,-10) h' 1 Mpc. It is notable that this 
is a very rough approximation of the true matter distribution - 
the IRAS map in Sigad et al. is heavily smoothed (12/i _1 Mpc 



Gaussian filter), and the density profile Jiras ( r ) ma t we use is 
measured by eye from their figure. However, it serves well for 
the accuracy that we aim for with the Tolman-Bondi method. 
The density models and corresponding cumulative masses are 
shown in Fig. |[ 

Figure [R^ shows schematically the system under study. A 
galaxy is located at distance d from us. The lines of sight to- 
wards PP and the galaxy form an angle 9. The distance from 
PP to the galaxy is then 



R 2 



d 2 - 2R PP dcos8. 



(12) 



The observable (radial) component of the peculiar velocity, 
^ci.obs (Eq. 0), of the galaxy is connected to the PP velocity 
field v r (r) by 



Vd.ohs = -v r (R PP ) cos 6* ± v r (r)] 



<,p P . 



(13) 



where + is used for d < Rpp cos 8 and — for d > Rpp cos 8. If 
we have a priori knowledge of the LG infall towards PP, v m f = 
v r (Rpp) we can solve PP peculiar velocities from observable 
quantities: 



iy(r) = ± 



Vd.obs + Vini COS 8 



1 



R 2 p sin 2 1 r 2 



(14) 
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Fig. 9. Assumed spherically symmetric density distributions 
<*>i (r), 62 (r) and Jiras ( r ), an d the corresponding masses inside 
radius r from the PP center. Overdensity 81 (r) and mass Mi (r) 
are plotted with continuous, (^(r) and M2(r) with dashed, and 
(5iRAs( r ) an d Mihas(^) w i m dashed-dotted lines. See text for 
details. 
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Fig. 11. LEDA galaxies within angular distance of 40° from the 
PP center. LG is at the apex of the cone, the gray sphere indi- 
cates the central PP region. The figure does not give a truthful 
representation of the galaxy distribution: Distances based on 
redshift velocities are distorted, selection effects are strong in 
larger distances and projection effects in the regions close to 
the cone axis. The figure is a projection perpendicular to the 
Galactic plane, which is marked with a dashed line. 
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Fig. 10. Galaxy distances from the PP and LG, and observed 
and PP rest frame velocities are plotted in this schematic pic- 
ture. Distances are normalized to the PP distance; Rpp = 1, 
i? s w 0.3. 



This is the equation used for a comparison of the TB models, 
giving the left hand side, and KLUN data, providing the right 
hand side of the equation. 



5. Comparison of Tolman-Bondi and KLUN sample 
peculiar velocities 



5.7. KLUN-PP sample 

Peculiar velocities v r of the KLUN galaxies are obtained in the 
following manner: 

- galaxies are binned into V intervals and distances d of the 
bins are calculated 

- each individual galaxy is given the distance of its bin 

- distances from the PP, 7-, and peculiar velocities, v r , are 
determined for each galaxy 

- galaxies are binned into r intervals and velocities v r of each 
bin are calculated 

First we extract a subsample (KLUN-PP) of KLUN galax- 
ies having angular distances from the PP center less than 40° 



(Fig. 11). For these galaxies we calculate distances d using the 
NDM and the mean surface brightness dependent diameter TF 
relation (Sect. 22). Strictly speaking, the bias correction pro- 
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cedure provides us with an unbiased average distance at a given 
d n or, as it can be shown, at a fixed kinematical distance dkin 
or corrected radial velocity. Thus we bin the galaxies into ve- 
locity intervals of AV^ = 500 km s _1 , and assign each galaxy 
the distance of its bin. 

To avoid any difficulties in combining possibly incompat- 
ible individual distance measurements we do not use ordinary 
weighted averages in the bin distance determinations. The dis- 
tances are deduced with more sophisticated Bayesian meth- 



ods (Press 1997 1: Probability distributions of "good" and "bad" 



measurements are calculated as 



and 



1 



27T(7logd,i 



■ exp 



(logdf - log ri g 

2 CT logd,i 



1 f (log^-logc^) 2 

exp 1 



2tt5 



2S 2 



(15) 



(16) 



where d' B is the (trial) bin distance, and di are the individual 
galaxy TF distances with errors cri og d,i- These individual errors 
are calculated by 



'TF 



'log D, 



(17) 



where cttf = 0.125 (Theureau 1998) is the TF relation scatter. 



and ciog £> 25 i are the individual (log) apparent diameter mea- 
surement errors. In Eq. (|l6|) S characterizes the standard de- 
viation of "wrong" measurements. In practice S is assigned a 
large enough value, exceeding the range where the measure- 
ment error would go unnoticed. We use here S = 0.2 (in units 
of log d). 

The probability that the binned set of data (D) have a dis- 
tance d' B is given by 

P(logd B |X>) oc /j] [pP gii + (1 -p)P h ,] dp (18) 



following Eq. (16) in Press (1997). Here pis the probability that 
individual distance is calculated correctly (and thus 1 — p is the 
probability for a wrong distance determination). We sum over 
all combinations of individual distance determinations being 
carried out correctly/incorrectly and over all values of p. As- 
suming uniform prior probability distributions for p and log d^, 
the sum reduces to Eq. (jilt). The proportionality in Eq. ( |l8| ) 
hides the constant needed for Yliogd' PQ°gd B \V) = 1- 

Actual bin distance c2b is adopted from the peak of the 
P (log d' B \T>) vs. logd' B curve, and 95 % confidence levels 
are deduced from calculating the area below the curve (Fig. 
|l2| ). Figure [l3| shows the individual and binned values for all 
the data points in the KLUN-PP sample. For comparison the 
binned distances are calculated also using the simple weighted 
mean, see Fig. |l2| There are no significant differences to the 
values obtained with the bayesian method, the two methods de- 
viate < 1 a for each d B - 

The peculiar velocity field around PP is determined simi- 
larly. Firstly each galaxy is assigned its bin distance (de) and 
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Fig. 12. An example of using the P(log d B |2?) vs. d' B curve in 
the bin distance determination. Individual galaxies are plotted 
on the d' B axis. Vertical lines show the adopted de and the 95 % 
confidence levels. For comparison, the simple weighted mean 
bin distance is also drawn with 2 a error bars. KLUN-PP galax- 
ies with 2500 km s _1 < V < 3000 km s _1 are shown in this 
example. 
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Fig. 13. KLUN-PP galaxies, bin distances and 95 % error bars 
in a V vs. d diagram. 

95 % limit errors (denoted by ob)- PP distances (ri) are calcu- 
lated by Eq. (|l2|), and errors by 



max(|rj - r(d B + er B )|, \n - r(d B - (Jb)|) 



(19) 



where r(x) is the function of d given by Eq. (|l2|). Thus Eq. 
( |T9| ) gives a safe estimate for the maximal error in the 95 % 
confidence interval, unless the galaxy has 



<7b < Rpp cos 9 < d B + &B ■ 



(20) 
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Fig. 14. Binned KLUN-PP data points and TB v r (r) curves. 
KLUN points are fixed with v- in { = 0. The error bars are 2 a 
confidence intervals. TB curves are calculated for the three 
density models, Si, &2, and 5tras- For eacn density, there is a 
set of curves, corresponding to the seven different cosmologi- 
cal models, (fi , Q A ) = (0.1,0.9), (0.1,0.0), (0.2,0.8), (0.2,0.0), 
(0.4,0.6), (0.4,0.0), and (1.0,0.0). 
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Fig. 15. The formal 1 a confidence regions (see text) for vi n f 
and Mpp. The best fit values are marked with triangles for 
= 1 — Oo and circles for fl\ = 0. The dotted lines show 
the Mpp limits obtained in Sect. H. 



The point at d = Rpp cos 6 has the lines of sight towards us and 
PP exactly perpendicular. There \r% — r(d)\ has a local maxi- 
mum. At this point the peculiar velocity towards PP is also hard 
to estimate (a v goes to infinity!). Thus we reject all galaxies re- 
stricted by Eq. (Eob. 



Peculiar velocities are obtained with Eq. (14) and errors by 



a Vt i = max{|t>i - v r (x)\;x G (d B - cr B , d B + cr B )} (21) 

giving again a safe 95 % confidence level estimate. In the end 
there are 220 galaxies with PP-distances and peculiar velocities 
from 354 galaxies in KLUN-PP before the restriction of Eq. 
(20). From these, 180 have r > R s . Thus roughly half of the 
original data are lost due to the constraints. 

Galaxies are then binned into PP distance intervals of Ar = 
0.05 Rpp. For each bin we calculate PP distances and peculiar 
velocities with 95% confidence levels using the Bayesian meth- 
ods stated above. Notice that for the peculiar velocity calcula- 
tions we need to assume a value for our infall velocity towards 
PP. That is, we have v- m { = v r (Rpp) in Eq. (p"4|), as a free 
parameter. 

5.2. KLUN-PP vs. TB 

We study how well our TB models produce the velocity field 
of the KLUN galaxies when the models are varied by 

- f2o, for which we consider values 0.1, 0.2, 0.4, or 1. Major- 
ity of recent studies agrees with this range of values (Pri- 



mack 200C) 



- a choice of either a flat A ^ (i.e. f^A 
open A = universe 



1 — f2rj) or an 



- a choice of density distribution around PP, assumed to be 
one of the toy models, 5\ or 62 (Eqs. |l0|-|lT|), or the spher- 
ically symmetric approximation of the IRAS density map, 

<5iras 

- mass of the central part of PP, for which we assign values 
between 2-15 /i _1 10 15 M Q . In the case of ^iras th e mass 
is varied by a linear bias factor. 

KLUN-PP data are adjusted by 

- LG peculiar velocity towards PP, v- m f (Eq. 113), for which 
we give values between and 600 km s~\These limits 
were chosen because km s _1 is a natural lower limit for 
a TB application, and after all we expect a large mass con- 
centration such as the PP supercluster to have at least some 
influence on the LG. In the other end, velocities exceeding 
600 km s _1 would have already been noticed in our motion 
with respect to the CMB, and are therefore not studied. 

Figure [14] shows the comparison of the peculiar velocities 
observed for KLUN-PP data and those predicted by the TB 
models. Binned KLUN-PP data points with 95 % confidence 
level error bars are plotted in the case of u; n f = 0. The curves 
are best fit TB curves to KLUN points, obtained by varying 
Mpp in the TB calculations. Two things are immediately ob- 
vious from looking at Fig. [IJ: Firstly, the IRAS density model 
fits much better to the observations than the toy models Sec- 
ondly, different cosmological models give very similar best fit 
curves. 

To put this into numbers, we describe the goodness of fit by 



N 

x 2 = ^ 



i=l 
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Table 4. The best fit values and the formal 1 a confidence re- 
gions of Mpp and &iras- See text for details. The values of 
Mpp are given in units of h~ x 10 15 M c . 



' ■ ■ 



0.1 

0.1 
0.2 
0.2 
0.4 
0.4 
1.0 



0.9 
0.0 
0.8 
0.0 
0.6 
0.0 
0.0 



2.7 
3.1 
3.9 
3.9 
5.5 
5.5 
9.4 



Mpp 

TTT 



0.6 
+1.4 
1.0 
+1.8 
1.4 
+ 1.8 
1.4 
+2.2 
1.4 
+2.2 
1.4 
+2.9 
1.8 



0.47 
0.41 

0.49 



frlRAS 
+0.U8 



-0.14 
1+0.14 
-0.11 
1+0.24 

-0.15 

49+ 24 

"•^-0.15 
j+0.20 
-0.16 
j+0.20 
-0.16 

u - dd -0.17 



0.53Z 
0.53;! 



where Vi is a KLUN data point, Ui its 1 a deviation, approx- 
imated to be half of the 95 % error bar, and Vi is a point 
from a TB curve, at the same r as the KLUN point. We get 
X 2 = 8.67 ± 0.20, 12.56 ± 0.37, and 2.01 ± 0.11, and sig- 
nificance levels P = 0.796 ± 0.015, 0.482 ± 0.029, and 
0.99976 ±0.00007 for 5 U S 2 , and S IRAS , respectively. The lim- 
its reflect the variance between different cosmologies. 

Figure |l5| shows the formal 1 a confidence regions for the 
values of Vi n { and Mpp obtained from fitting the (Jiras TB 
curves to the KLUN data points. The confidence regions are 
defined in the following manner: The value of \ 2 is calculated 
for each pair of parameters (vi n f, Mpp). Then a point in the 
(v,M) space belongs to the confidence region if x 2 (v,M) < 
y 2 ; „ where y 2 . is the minimum value of y 2 in the whole 
region of v = 0-600 km s" 1 and M = 2-15 /i" 1 1O 15 M , and 
a v is the average of the 1 a errors of the binned KLUN ve- 
locities. This definition thus gives all the points in (fi n f,Mpp) 
giving x 2 ~ Xmin within observational errors. 

The confidence regions in Fig. [l^ suggest that the LG infall 
towards PP is less than about 100 km s . In each case the best 
fit, plotted in the figure as a black triangle for Q,\ = 1 — fio 
and a white circle for fiA = 0, is obtained with t>i n f = 0. 
The best fit values of Mpp are coupled to the value of fio; 
the correspondence with the range Mpp = 4-7 h^ 1 10 15 M©, 
which was obtained in Sect. is better for fio = 0.2-0.4 than 
for fi = 0.1 or 1. 

As mentioned above, the mass Mpp is varied by adjusting 



a bias parameter for the <5-values given by Sigad et al. (1998). 
Let 5s be the overdensity at different distances around the PP 
center, obtained using the contour plot in Fig. 6 of Sigad et 
al. (1998). Since Sigad et al. use fio = 1 in their figure, <5s(r) 
must be multiplied by fi^ 0,6 for a general case of fio 7^ l.Then 
we assume a bias of the form fip ' 6 ^ = ^iras £>iras> where 
(Jiras gives the density which we use in the TB calculations. 
Increasing &iras would mean decreasing the value of Mpp, 
for example when fio = 0.4, Mpp = 4 h~ 1 10 15 M© requires 
friRAS = 0.83 and M PP = 1 hr 1 1O 15 M means 6 mA s = 
0.42. The best fit values for Mpp and 5iras> an d their formal 
1 a confidence regions, derived from the regions in Fig. [fjl are 
listed in Table (§ 

Our conclusion from comparing the KLUN-PP peculiar ve- 
locity data to the TB-modeled velocity field is that we get a 
good correspondence between the two when 
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Fig. 16. Binned data points in v r vs. r diagram after the itera- 
tive process for correcting distance determinations. Each time 
the PP distances are calculated, approximately half of the data 
are lost due to restrictions (Eq. |o|) and r > 0.3 i?pp. Only six 
binned data points are now present, rather than the 13 shown 
in Fig. |l4| In this figure we have adjusted the KLUN points 
with v- m i = 0, and the curves show the best fitted TB solutions 
using ^iras- The best fit values are Mpp = 2.9, 3.7, 5.5, and 
9.4 xh- 1 1O 15 M , for fi = 0.1, 0.2, 0.4, and 1.0, respec- 
tively, corresponding to bias values 6iras = 0.45, 0.52, 0.53, 
and 0.55. 



- the density field around PP is approximated from the IRAS 
galaxies, as presented by Sigad et al. (1998), with addition 
of a linear bias &iras ~ 0.5 

- the infall velocity of the LG towards the center of PP, t>i n f , 
is less than 100 km s . 



We also note that 

- varying the value of fio between 0. 1 and 1 and fixing fiA = 
1 — fio or have no significant effect on the results, except 
that 

- the mass of the central part of the PP agrees with the limits 
obtained with the virial theorem if fio = 0.2-0.4. 



5.3. Iterative corrections, p-class analysis 

The first results can be used to correct the observed velocities 
in our NDM distance determination. With the best fit <5iras TB 
model shown in Fig. |l4] we calculate peculiar velocities at dif- 
ferent distances r, which are used to correct the input kinemat- 
ical distance scale. The whole process of binning and solving 
for PP distances and peculiar velocities is then repeated with 
these new kinematical distances as a reference scale. 

A big disadvantage in this iterative procedure is that at each 
iteration we lose about half of our data points due to the con- 
straints of Eq. ( ^o[ ) and r > R s = 0.3 i?pp. Figure |l6| shows 
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how we have only 6 binned data points containing 90 individ- 
ual galaxies after the first iteration. Our conclusion is that the 
iterative process does not in this case improve the results. 



and looking at the galaxy distribution in Fig. |[ PP hardly can 
be taken as such. However, a smoothed density map, like the 
IRAS map by Sigad et al. (1998), has a spherical structure. 



Ekholm & Teerikorpi (1994) used in their study of the pe- 
culiar velocities in the Great Attractor region a so-called p- 
class analysis. There the possibility of the log c?tf not being 
normally distributed is taken into account by taking only the 
galaxies lying in the unbiased plateau of the (\ogHo)(d) dia- 
gram. Applying the same procedure here did not give disagree- 
ing results with the velocity patterns shown e.g. in Fig. |l4|, but, 
again, the number of points given by the p-class analysis was 
small. 



6. Summary 

We have studied the velocity field around the Perseus-Pisces 
supercluster using the KLUN galaxy sample and Tolman- 
Bondi modeling. In the TB calculations we used a sm oothe d 
density field of IRAS galaxies, as given by Sigad et al. ( 1998 ). 
There is a good correspondence with the peculiar velocities 
given by the KLUN data and the TB models, if the infall ve- 
locity of the LG towards the PP center is assumed to be less 
than 100 km s _1 and if we add a bias in the Sigad et al. density 
field, fciRAS ~ 0.5. The resulting small infall velocity of the LG 
towards the PP region, t>; n f < 100 km s _1 disagrees with the 
v i n f « 350 km s _1 of the earlier kinematical studies (see Sect. 
if), bu t is consistent with the recent POTENT analysis (Dekel 
1999|). 



The goodness of fit between the KLUN and TB velocities 
is not affected by varying the cosmological parameters flo and 
fl\ within reasonable limits. However, requiring the mass of 
the central region of PP to be in the limits obtained with virial 
calculations, Afpp = 4-7 hr 1 10 15 M Q , we get constraints for 
the value of VLq. Then the intermediate values £7o = 0.2-0.4 are 
preferred to the more extreme fi u = 0.1 or fin = !■ The value 
of Ha, assumed to be either 1 — or 0, does not have a sig- 
nificant effect on the results. In fact, the infall velocities around 
growing structures being little influenced by the presence of A 



is not so unexpected (Peebles 1984, Lahav et al. 1991 ) 



Calculation of the virial masses of the most prominent 
clusters in the PP region resulted in cluster mass to lu- 
minosity ratios M/L = 200 - 600 h Mq/Lq. Combin- 
ing this with a luminosity density of the universe L = 
10 i9.6±o.i h w Hz -i Mpc -3 (Li ny et al. |199§), and assum- 
ing that the M/L of the PP region is valid in the rest of the 
universe, we obtain = 0.1 - 0.3. Allowing M/L to grow at 
larger scales would mean that this is a lower limit. The result 
agrees with recent estimates of f2o- In comparison, in the study 
of the Hubble flow around the Virgo cluster, Teerikorpi et al. 
( 1992 ) found that the constraint M/L = constant is inconsis- 
tent with Slo < 0.2 or f2o > 1.4. A smaller Slo was possible if 
M/ L in Virgo is larger than in the field, contrary to the usual 
bias scenario. The same is valid here. 

One of our aims was to study the applicability of the TB 
method in a complex environment like the PP supercluster. TB 
solutions are for spherically symmetric mass distributions only, 



Similarly, the small scale dynamics are smoothed away, when 
we use an averaged velocity map, where the radial parts of the 
peculiar velocities are averaged over a shell around the main 
galaxy concentration. This is shown not only in the compatibil- 
ity of the TB results and the KLUN data, but also in a study of 
an N-body simulation. There it was confirmed that the radially 
averaged velocity field can be calculated by the TB method, us- 
ing the observed smoothed density field. There were marginal 
differences between the TB solutions and the simulation when 
we looked at the peculiar velocities around large clusters at dif- 
ferent epochs. However, at the present epoch these deviations 
are negligible. We argue that the Tolman-Bondi method can be 
used for a reliable deduction of a radially averaged velocity 
field around large galaxy concentrations. 
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